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THE  REPRESENTATION  OP  AN  EMPIRICAL  FUNCTION 
BY  MEANS  OP  A POLYNOMIAL. 


by 

D.  J.  Behrens 


Abstract 


If  the  values  of  a function  are  given  at  n equally- 
'spaced  intervals  of  the  argument,  a polynomial  of  order 
n — 1 can  be  found  which  corresponds  exactly  to  the  given 
values.  In  general  a polynomial  of  order  r(  < n-1) 
cannot  be  made  to  do  this,  and  ;Jj.('f),  the  minimized  sum 
of  the  squares  of  the  deviations,  will  exceed  0. 

It  Is  shown  in  this  paper  how  j|j,(  f),  and  the  leading 
coefficient  of  the  polynomial  in  question,  can  be  evalu- 
ated without  having  to  solve  any  'simultaneous  equations. 
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1.  introdaction 


It  frequ.ently  happens  in  experimented  work,'  that  we  can  measure  the 
value  of  some  fupction  f at  several  values  of  the  argument  x, ■ and  we  wish 
to  find  a polynomial  p(x)  vhich  is  a sufficiently  accurate  measure  of  the 
given  functional  values.. 

We  shall  confine  our  attention  to  the  case  in  which  (1)  equ,al  errors 
are  considered  to  be  equally  probable  in  each  of  the  given  functional 
vedues,  ■ so  that  the  method  of  least  squares  may  be  ujsed  without  weighting 
factors  to  give  the  "best"  polynomial,  and  in  which  (2)'  the  values  of  the 
function  are  given  for  equ.ally  spaced  values  of  the  argument  x.  We  shall 
assume  them  to  be  given  at  x = 1, 2,  k...n,  and  the  corresponding  functional 
valu.es  will  be  denoted  by  fi  fj  .....  f». 

The  two  reservations  made  in  the  previous  paragraph  can’  upu.ally  be 
ensu.red  by  suitable  experimental  techni'qu,e.  We  make  them,  because  on  their 
validity  depends  the  simple  method  developed  below. 


It  is  clear  that,  if  n functional  valu.es  are  given,  we  can  construct  a 
polynomial  of  order  n-1  which  will  pass  exactly  through  each  point. 

In  general,  a polynomial  of  order  r,  less  than  n-1, -will  not  pass  exactly 
through  all  the  points.  We  can  however  define  an  r^^  order  polynomial 
Pj,(x)  by  adjusting  its  coefficients  so  as  to  minimize  the  sum  of  the  squares 


^ f 1 ® 

This  minimized  sum  will  be  denoted  by  J*)j,(f). 

I.f  now  we  examine  the  fom  of  as  r is  increased  from  -l*to  n-1,' 

we  find  that,  if  the  functional  valu.es  ft  f2. , . .f„  really  do  lie  close  to  a 
polynomial  of  comparatively  low  order,  the  general  form  of  ij.(f)  is  as 
sketched  in  Pig.  1. 


In  the  region  shown  as  (a), -we  have  too  few  parameters  to  represent 
adequately  even  the  general  ou.tlihe  of  the  cu.rve,  ■ while  in  (o)  we  have  so 
many  that  phoney’  "accuracy"  is  obtained  by  fitting  a cu.rve  even  through  the 
random  errors  of  observation.  Ihe  region  (b)  is  that  in  viiich  we  have 
enough  parameters  to  represent  the  general  form  of  the  cu.rve,  bu.t  not  so 
many  as  to  purport  spu.ri6u,s  accu.racy.  The  valu.e  of  r in  which  we  are 
interested  lies  at  the  lower  end  of  this  region’,  and  is  indicated  by  the 
dotted  line  in  Fig.  ;1. , 


We  shall  show  how  to  find  thi's  vadu.e  of  r withou.t  the'tedioup  and 
heavy  arithmetical  work  involved  in  actually  calculating  p{x)  for  each  r. ; 


Fi'gu.re  1. 


log  <Sj. 


y « constant  i s a polynomi al  of  order  zero:  - i t i s conveni  ent  to  define  y « 0,  efiich 
has  one  fener  degrees  of  freedom,  as  a polynomial  of  order -1. 


•1* 


8.,  Determination  of 


Let  Aj  denote  ~ fji'  Aj  - Aj - Aj , and  in  general  aJ  « A^-Aj— 
TMa  is  the  usual  fihite  -difference  notation'  for  forward  differaices. 

Let  Itlj.  denote  the  mean  of  the  r^  order  differences  A*^,  weighted  in 
accordance  with  the  formula 


where  the  terms  in  roufid  brackets  ( ) daiote  binomial  coefficients.  [It 

will  be^proved  in  appoidiic  (proof  1)  that  "F  in  fact  equal 

^ (gr+i)  ’ ^r  defihed  as  above  is\rr  fact  a Weighted  mean  of  the 

a'^'s.;] 


Itlo  is  defined  ih  the  same  way,  ranonbering  that  Aj  = fj,  and  in  fact 
Itlo  can  be  seen  to  be  the  arithmetic  mean  of  the  f's.  In  general,  lUj.  is 

the  r''^  difference  of  Pj.(x),  the  r^^-order  polynomial  gi'vihgthe  closest 
fit  to  the  given  f's  (See  Appendix,  ■ Proof  8). 


Th'en  f)  and  Aj>-]_(  f)  are  connected  by  the  relation' 


Vi(f)  7 Vf) 


See  Appendix,  > 
proofs  3 and  4.; 


Sihce  we  know  that  4,-1  ' 0,-  and  that  i4_i(f)  » 2 (ff),'  we  can 

calculate  the  «j,'s  beginning  from  either  end,'  Alternatively' the  fact 
that  both  end-valu,es  are  known  serves' as  a check  on  the  ariithraetib  if  we 
calculate  f)  - 4^ -)  every  r from  0 to  n-1. ; 


3.  Probable  error  in  til j,. 

It  ronaihs  to  calculate  the  probable  error  of  the  Hi's,'ih  terms  of 
that  of  the  giver  functional  valu.es. 

Erqrreasihg  tHj,  ih  terms  of  the  f's  themselves,  instead  of  their  r^^ 
differences,  ■ we  get 


r 

2 

m-o 


(-U*  C) 


P) 


And  thup  the  mean'  squ.are  error  in  lUj.  is  given  ih  terms  of  that  ih' 
the  fis  byi 


6“  lllr/6»f 


n 

2 


Li  o 

(**i 

fn-j-m+r^l 

\ 

f n+r  N 
^2r  + lJ 

.2. 


This  redupes  to 


6“  tHj,  / 6®f  = 


[see  Appendix,  proof  s] 


Bu,t  this  ratio  has  alreac^  been  shown'  to  be  equ^  to  /(Ajwj^(  f)  - Ajj(  f) ), 
so  we  have  the  result  that 


6®  Itl^  6»f 


Often,  we  do  not  know  the  mean  squ,are  error  6“*^  We  may  however 
assume  that,  liif  r has  been  sensibly  chosen,  it  is  of  the  order  of  i»p(  f)/(n-r-l) 
- we  divide  by  (n'-i“l)  as  this  is  the  excess  number  of  degrees  of  freedom  over 
those  recjii'red  to  fix  the  parameters  of  Pj-Cx). 


Thup  we  may  say  = )ij.(  f)/(n-r-l)()lij^2(  " V 

condition  that  should  be  significantly  different  from  zero  is 
A 'f)  < (n-r-l){iJj^l(f)  - \(f)>, 


or 


f ) 


(n-r-1) 


f ) 


(n-r) 


It  is  therefore  worth  while, 'in  tabuj.atihg  Aj-lf),'  also  to  note  the 
valu.e  of  the  quptient  i)j.(  f) /(n-r-1),  since  (a)  when  6*f  is  known,  it  is  this 
qu.antity  which  must  be  kept  dowi  to  the  order  of.6®f,  and  .(b).  in'  any  case  if 
the  quptient  ^ip(  f )/(n-r-l)  exceeds  the  corresponding  valu.e  at  (r-1),  this 
shows  that  no  significant  additional  accu,racy  arises  from  taking  an'  r''^  order 
polynomial  instead  of  Pj-p;  Ibis  argument  is  used  in  Example  3 below. 


Reference  to  other  work 

M-  G.  Snith  "The  Lorentz  Method  of  Analysis  of  Experimental,  Data 

using  Orthogonal  Polynomial's"  (A-R-E-  Report  No.  15/58). 

[This  paper  presents  a method  for  calculating  the  best- 
fitting r^^  order  polynomial,  having  first  obtained  that 
of  order  r-1.  It  differs  from  the  treatment  of  the 
present  paper  in  computing  first  the  polynomials  than- 
selves,  and  only  afterwards  the  values  of  ii(f),  whereas 
we  have  here  treated  the  residual  errors  as  being  of  prime 
interest  - i.e.  we  determine  the  required  order  of  polynomial 
without  having  first  computed  it.  Either  treatment  can 
be  more  useful  than  the  other,  according  to  circumstances.  ] 


EXAMPLE  I . Functional  valu.es  lying  close  to  a polynomial  of  low 

order  (cu.bic). 


£61  3;'8  £02  320  198  63  -69  -162  -181  -100 

1111111  1 1 1 
£61  378  392  320  198  63  -69  -162  -181  -100 

~Ii7  14  ^2  -135  -132  -93  ^19  E 

9 16  21  24  25  24  21  16  9 

1053  224  -1512  -2928  -3£75  -3168  -1953  -304  729 


-103  -86  -50  -13  3 39  74  100 

86  84  126  150  150  126  84  36 

-3708  -7224  -6300  -1950  450  4914  6216  8600 


17  36  37  16  36  35  26 

84  224  350  400  350  224  84 

1428  8064  12950  6400  12600  7840  2134 


19  1 -21  20-1-9 

128  350  525  525  350  126 

2394  350  -11025  10500  -350  -1134 

Ea  ^ 41  ^^1  ^8 

ise,  336  441  336  126  " 

-2268  -7392  13081  -7056  -1008 


-4  63  -62  13 

84  196  196  84 

-336  12348  -12152  1092 


67  -125  75 

36  64  36 

2412  -8000  2700 


-192  200 

9 9 

-1728  1300 


165 

11234 


1716 

51466 


136 

-2388 


= 584008 

4 = 584006-  (1100)^(10x1)  = 5840(»-121000  = 463006 
Si  = 46300e-(H234)7(ie5x2)  = 463OO0-3824S2®7ier8CB75®7i 


58401 
? 4 ~ 51445 
■5  ;4l  ~ 1QJ?2 


Si  = 805?5"7io5-(4D02)7('?92x6)  = 8O0?5°7iof3STO‘'7i3s=778O5*78no  ^ ^ " 11C89 

S,  = T?205*766o-(5i4K)7(17l6xa0)«7;2C6“7eeo-77178**/e»Bo=27“'/Bse  - 4.5 

4'  = 27*/£,sB-(9S6)7(S0C2)W)«87“79Bo-3*°7i72  » 23*7io2  i 4 ~ 4.6 

4'=  23*7i3i-(357)7(ia65«58)=23“7i»a-*"/T8o  * 22"**/ai*5  ^4-5.7 

4 - 22”7ai46-(958)7(50O)«01)  * 28”*72i*5-”*/7eo  » 2l“7*i45  i 4 ^ 7.'0 


2l”7«48-(288e)7(136>0432)-El”7«*«-17j^^  = 3”y«as 
3‘’V9  3S-(72)7(18x32er70)»3”V.«r"77«  = 3“*Vi2«b 


■i  4 ~ 7.'0 
i 4 ~ 1.6 


3*“Visi6«r(3ffi)7(lx48K0)  = 3“*  Vi«8«r3“®  Vuibb  = 0 


A 


EXAMPLE  II 


Functional  Values  Not  Lying '.Close  to..a  Polynomial  of  Lew  Order. 


Independent  variable:- 

0 

1 

2 

3 

4 

5 

6 

SUMS 

functional  values  :~ 

1 

2 

4 

8 

16 

32 

84 

weighting  factors 

1 

1 

1 

1 

1 

1 

1 

7 . 

ttlo=  7 

1 

2 

4 

8 

16 

32 

64 

127 

first  differences 

1 

2 

4 

8 

16 

32 

weighting  factors 

6 

10 

12 

12 

10 

6 

(IImS  n n 

6 

20 

48 

96 

160 

192 

522^ 

second  differences 

1 

2 

4 

8 

16 

weighting  factors 

15 

30 

36 

30 

15 

126y„  ilE. 

1^2“  42 

699 

15 

60 

144 

240 

240 

third  differences 

1 

2 

4 

8 

weighting  factors 

20 

40 

40 

20 

120.  z 

lUo*  2 

20 

80 

160 

. 160 

420 

fourth  differences 

1 

2 

4 

weighting  factors 

15 

25 

15 

nl4=  1 1 

15 

50 

60 

125 

fifth  differences 

1 

2 

weighting  factors 

6 

6 

12».  1 

6 

12 

18  ' 

sixth  differences 

1 

weighting  factors 

1 

^ Ule=  1 

i>»  = 0 


iig  ~ lie 
li*  ~ lie 

— i>4 
is  — i>3 
ii  — is 
io 

i-i~  io 
Check:- 


II 

0 

92  4 “ 

1 

824 

1 

^5  ® 924 

1.1  X 10"" 

12 

S«2  * 

12  '9 

2S2  4 ■ 

9 

2 e 

13  1*00  25 

®4-  92  4+2  8=  92  4® 2 3 1 

26 

^ = 231 

'K/ 

1.1  X 10"‘ 

^ V.2 

70  Xffl'4  “ 

00  0 8 0 
7C  ^12 1 ” 

02  5 

' 154' 

25  025  192  5 20 

®3  = 231+16*=  402  = 0' 

^ = T 

'^w 

4.2 

1^0  /7\2 

20^\2  J 

147 
= 2 

2 5 147  406_ 

*^a=S+2=0  ® 

293 

^ = “5“ 

<"0 

7.8  X 10 

12Q  /233 

0 (v  42 

}"=(233) 

Vs4 

„ . 239^  20271 

i>i=233{^  + 7r}=TF 

202  71 
^ - 26 

7.2  X 10= 

Be  /®0i\ 
2 \ 2 8 / 

* =(261) 

00 

/4o«  (2CE71+68]21)/28=220G^7i 

220  98 

ib  = 7 “ 

3.2  X 10® 

f (x) 

“ =(127) 

V7 

(I6129*-2209e)/7=5461 

iLi=  5461 

JV 

5.4  X 10® 

Sum  of 

f“  = 1 

+ 4 + 

16  + 64  + 256  + 1024  + 4096  = 5461 

10 
1 

0.1 

And  the  shape  or  the  curves  0.01 
suggest  that  the  functional  values 


and  V(n-r-iX 


certainly  do  not  fit  any  loiv-order 
polynoml  el. 


0.001 


'X  , 

"o-.. 


m|  u 


EXAMPLE  III 


What  is  the  lowest  order  of  polynomial  consistent  with 
the  following  observations, 


(a) 

if  the  probable  error  of  each  observation  is 

50 

(b) 

II  II 

II  II  II 

It  II  II 

1 

r X 1 2 

3 

4 5 6 

7 8 

9 

Lf(x)  27  676 

1211 

1576  1730  1640  1309  775 . 

116 

f 

27  676  1211  1576 

1780 

1640  1 309  775  116 

E(A  wf) 

wf 

1111 

1 

1 1 

L 1 

9060 

L 

649  535  S65 

154  -90 

-331  -534 

-659 

8 14  18 

20 

20 

18  14 

8 

1186 

-114  -170  -211 

-244 

-241  -203  - 

-125 

wf 

28  63  90 

100 

90  63 

28 

- 95271 

A® 

-56  -41 

-33 

3 

38  78 

wf 

56  140  : 

300 

200  140  56 

- 5188 

A* 

15  8 

36 

^ 40 

wf 

70  175 

225 

175  70 

19475 

A® 

-7 

28 

-1 

5 

wf 

56  126  126 

56 

3290 

A* 

35 

-89 

6 

wf 

28 

28 

- 873 

A^ 

—64 

35 

wf 

8 

8 

- .232 

A® 

-99 

wf 

1 

- 99 

is 

= 0 

is  = 0 

i,  = (99  ) 71  X 12870  = 

0.76« 

is  = 0.76 

i,  = 0.76 

= (.87)* 

is 

- i,  = (232)716x3432  = 

0.98 

is  = 1.74 

■■i  is  = 0.87 

= (.93)* 

is 

- is  = (273)7105x924  = 

0,77 

is  = 2.  51 

is  = 0.84 

= (.92)* 

■i4 

-is  = (3290)7364x252  = 

118 

i*  = 121 

is  = 30 

= (5.5)* 

is 

- i4  = (19475)7715x70  = 

7578 

is  = 77  X 10* 

-1  is  = 15. 4 X 10  * 

= 30* 

i* 

- is  = (5188)7792  x20  = 

1701 

is  = 94  X 10* 

■i  is  = 15.3  X 10* 

= 39* 

ii 

-is  = (95271)7462  x6  = 

32.7  xlO® 

i,  = 32.8  X 10® 

4 i,  = 4.7  X 10® 

= (690)  * 

io 

- ii  - (1186)7120  X 2 • 

5861 

io  ■ 32.9  X 10® 

i io  ■=  4.1  X 10® 

• (640)* 

- - (9060)  79  X 1 . 

91.2  xio" 

•i_j«  124.1X10® 

i i_^»  13.8  xio® 

• (1200)* 

.7. 


/Aj./(n-r-l)  shouJ.d  be  of  the  same  order  Of  magnitude  as  the  probable  mean 
error  of  observation  [for  detailed  treatment  see  any  textbook  of  statistics,' 
under  the  heading.  "I*  distribution' of  gooAiess  ot-,  fit"].. 

In  the  example  given,'  it  will  be  seen  that  a parabola  (r»2)  is  consistent 
with  a probable  error  of  50, 'While  if  the  probable  error  is  only  upity, ' a 5^^ 
order  polynomial  is  required  to  give  an  adecjr.ate  fit.. 

With  a probable  error  of  50  in  each  f,>  however, ' the  lower  valu.es  of 
/>^j,/(n-r-l)  for  r ih  excess  of  3 .vould  appear  highly  suspicious,  and  vould  lead 
very  definitely'  to  suppose  that  the  errors  are  largely  of  a systematic  rather 
than  a random  natu,re. 


.8. 


APPENDIX 


Proof  1. 


Y (Y)("r)  ■ (^;i) 


[n  > r],; 


Let  6(n,  r)  denote 


Y (Y)Cr)  ■-(-.) 


Writing  the  binomial  coefficients  as  the  susns  of  coeffici'aits  of 
lower  order,  ■ we  may  say 

- Y (Yf) (i  * ”-^)(Y?)(i * Y)  * lYi) 

Iheteim  I)  vanishes  both  at  j“l  and  at  j=n-r;  we  may 

_2 

therefore  shorten  its  range  of  supimation  by 


one  at  each  end. 


so  that  6(n, 


.J  * i J • ^ 


— 6(n-l,r-l)  + 6(n-2,.r) 
r 


Now  6(n,0)  = (V)('^o^)-(l 


= n-n  = 0 


Again;  at  r = n-1,  j.  can  only  be  1,  and  6(n,n-l)  = (n-l)(n-l)  (ai-1 ) " ^ ^ ® 

At  r • n-2  ,j.  can  be  1 or  2, 'and  6(n,n-2)=^n-2)(n-2)'^(n-^(|i^2)"^^2)°^^‘“^^^^  ° 
Ihis  suffices  to  prove  by  induction  that  6(n,  r)  ■ 0.’ 


.(35' 


! general  form  of  connection 
■between  contigupu?  6's 
{iie.  (n,  r)  with 
(n-1,  i“l)  (n-2,  r). 


n=o  on  vrfii'di  6=o 


.9. 


Proof  2. 

That  E ^)/  (sr+l  ) the  difference  of  the 

r^^  order  polynomial  giving  the  closest  fit  to  the  points  fi  f* -fn 


Let  the  f's  be  cho'sen  'so  that  all  but  one  of  their  r''^  order  differences 
vanish,  and  the  other  Is  uni  tv.  Let  hf  =1.  Af  “ 0 for  k f Suph  a 
set  of  f's  is  given  by 
■for  k ^ j r. 

Let  Ppix)  = Sf.  x'"  + x^^ 


Is  unity.  Let  aJ  = 1.  ^ “ 0 for  k f j.  Su^ 

. fo  = f.  ...  = fj,rt  = 0 and  fk  = 


, . ajx  ao  . 
We  are  to  choose  the  a's  so  as  to  minimize 


j-n—l 


E (aj.x''  +...+  aix  ■t-ao)*' + 2 |aj,x^  + aj_j^ 

-1  x=j■^r  L 


x'^^+. . . ajx+ao 


-(tl 


If  we  use  a notation  in  which 
thi's,  and  'say  that  we  have  to  minimize 


) 


0 for  all  g < h , we  may  simplify 


" r 

'E  J a_  x^ 
x=l  I 


+ ap.^  x^l  + + aix  + ao  - j* 


We  obtain  a set  of  simultaneous  equations,  from  which  ap_^  ap_g  ....  a^  So 
can  be  eliminated,  leaving 


^r 

%r'  Sgr-l 

S2r“2  • 

...Sr 

^8r^i^2r^2 

’S21-3  • 

• • ^r-1 

^2r-2%r“3 

S21W.4  . 

• ■ ^r-2 

Sf.  Sp_;[_ 

^r-2  ’ 

* • • So 

n 

E 

1 

xr  ( 

1^1  / 

S21-1 

Sgiw2  •••• 

■Sr 

n 

E 

1 

x^^ll 

'x-j-i 

S2r-2 

%i^-3  •••• 

■Sr-1 

n 

E 

1 

xi-2 

M-1 
{ 1^1 

)'S2r-3 

^2r^4  *••• 

Si-2 

n 

\ 

E 

1 

[L 

Siw-i 

Sjwg  .... 

So 

where  ^ denotes  (1*^  ^ ^ n*^). 

We  have  to  calculate  the  values  of  these  determinants. 


10. 


1 1 may  be  shown  that 


1 

1 

1 

1 

= [3(r)]V3(2r+l) 

2r +1 

2r 

2r-l  "■ 

r+1 

1 

1 

1 

1 

2r 

2r-l 

2r-2  *" 

r 

1 

1 

1 

1 

2r-l 

2r-2 

2r-3 

rt 

1 

1 

1 

1 

r+1 

r 

r-1 

1 

where  3(r)  denotes  the  factorial  function  U 2.'  2.'  4.'  rJ 

From  this  it  follows  that 

= 3(n+r)  3(n-r-2)[3(r)J7[3(n-l)]*  3{ar*l) 


’S2r 

Sgr-i  • • • 

Sr 

S2r-1 

S2r-2' 

Sr 

■Sr 

Si-1  ••• 

So 

It  can  also  be  shown  that 


n 

2 

1 

Sgr-l  • • • 

■Sr 

= {..(7)(T)/(T)} 

-%r£'Sg[_Q 

Sp-l 

n 

y 

i^l/x-j-A 

Sgi^g  • • • 

Sr-l 

S2r-3  Sgr-4  * ’ ’ 

Si-2 

1 

^ V r-1  / 

Sp-l  Sp-g.i.. 

So 

n 

2 

1 

L 

CO 

So 

jy  3(n+r-l)  3(n-r-l)  3(r-l)  [3(r)]“ 


r 


[3(n-l)]  = 


3(2r) 


■go  that  a. 


(n-r-l)J  (2r+l)J 


(n+r)J 


And  hence  the  order  difference  of  PpCx),  v^ich  is  r.'aj,,  is  given  by 


C^r  ^)^(2r+l}  ’ fact  the  coefficient  of  aJ  i 


in  the 


expression  for  As  this  applies  for  any  j,  and  the  differences  and  Hi's 

are  all, linear  in  the  f's,  it  follows  that  Hip  is  always  the  r*^  difference 
of  the  best  pp(x),  whatever  be  the  values  of  the  f's. 


■i 


11. 


Proof  3. 


That 


f ) - f) 


is  a function  only  of  n and  r. 


Let  f denote  a set  of  fufictibnal  valu.es  fi  fa  fa  ....  f„-  given  at  n 
equally-spaced  iht'ervals'of  the  argutnent. 

Let  Pp(x)  denote  the  best  order  polynorai’al  through  the  functional 
values  - i'.  e.  the  parameters  are  so  chosen  that  t is  raihifni'zed  - and 
let  f ) be  u.sed  to  denote  this  minimized  su(n,  anA  f ) the  r^^ 

difference  of  the  polynomial  Pi-(x)., 

Now  let  f*  doiote  that  set  of  functional  falu.es  f*  f*  f*  suph  that 

the  first  (r+l)  of  them,  viz  f*  to  are  the  same  as  those  in  the  f set 

(f*  = fj,  f*  = fj,  ....f*.^2,)>  remaining  n-i— 1 values  are  so  chosen 

that  they  lie  on  the  r^^  order  polynomial  through  fi,  f*  ....;  fr+i* , 


Then  clearly f*)  •=  0. 

Next,  define  a new  set  f^0^  of  functional  values,'  su.ch  that  fj(9)  ■= 

9 fj+(l-9)f|, for  any  ,j  from  1 to  n.  It  is  clear  that  all  the  fj(0)'s  must 
be  linear  in  0,  • since  fj  end  f*j  are  constants  already  fixed.  From  this 
it  follows  that  ttlj,(f(0))  is  also  linear  ih  6,  ■ as  Itl-  is  linear  in  all  the 

Since  f(0))  is  linear  ih  6,  there  mu.st  exist  a 0 for  which  IUj.(f(0))  -0. 
Let  this  value  of  9 be  denoted  by  0**,  and  let  f**  represent  the  corresponding 
set  of  functional  values. 

Then  the  statement  that'ttVpl  f** ) = 0 means  that  the  best  r^*^  order 
polynomial  through  ft*' fj*  f^*  has  a leading  coefficimt  of  zero.  l.e. 
it  is  of  order  (r— 1)  or  lower.  Thu.s  it  followes  that  ^>j,(f**)  = i>jwj^(f**)« 

Now  the  coefficients  ih  Pj.(x)  [and  likewise  ih  Pjw]^(x)]  can  be  seen  to 
be  linear,  in,  the  functidncil  valu.es  and  hence  ih  9.  Thu.s  it  follows  that 
j^j,(f(0))  and  ^i~-i(f(e))  are  both  qu.adratic  ih  0.  The  difference 
iir_i(f(0))  - V^(0))  quadrati.c  ih  9.;  Bu.t  this  differaice  can  . 

never  be  negative  - (the  best  r^^  order  polynomial  mupt  always  be  at  least 
as  good  a fit  as, the  best  order).  Also  " '^r^f(9)^ 

vanishes  at  6 = 9**. 

...  Vl(^{0))  “ Vf(9))  = constant.- 

Again,  f0)  is  as  we  have  seen  lihear  ih  0,  • and  by  definition' vani'shes 

at  0 = 9** 

tkj,(f(0))  «=  K2(0  - 6**)  where  Ka  is  another  constant. 

From  this  it  follows-  at  once  that 

^r-1  ^^(0)^  ^ ^ ^ constant  independent  of  0.; 

tn;  h 

As  we  could  have  defihed  1=*  by- any  r+1  points  (not  necessarily  the  first  ones), 
we  can  see  that  this  constant’ mu. st  be  the  same  whatever  the  valu.es  of  the 
functional  valu.es, .and  ih  fact  r^f)  can' depend  only  on  the 

nupibers  n and  . , - - - - 


. 12. 


Proof  4.  Th.it  ^ is  equal  to 


NOTE:,  In  proof  (3)  above  we  have  shown  that 


- Vn 


tH*(r) 


depends  only  on  r and  a,,  and  not  on  the  f's.  Cbnsequently  we 
are  at  liberty  to  choose  now  what  f's  we  wish 


* * » 


Let  the  f's  be  chosen  ao  that 


fj  = for  1 ^ .j  i r+1, , 


and  fj  = 0 for  r+1  < j < n, 

n . 

Then  2 fj  = 0 for  k ? r and  = 1 for  k = r. 

■j-1  ^ 

n 

Prom  this  it  follows  that  Pjwj^(x)  » 0 and  f)  = 2 (fj 

Por  Pp(x)  we  have  the  conditions 

^r  ^2r  ^r"l  **"  ^r+1  "**  ^r  ” ^ 1 

i 

^r  ^2r-l  ^r-1  ^2 r-8 ^r  ^r-1  ' ^ 

^r  %r-2  ^r-1  ^2r“3  * a Sjw]_  + ao  Sj^-g  » 0 


'*^r  ^ ^r^l.  ^ ^ ^ So  So  * 0 J 

n 

• 2 (ap  + ap_j^  aij  + ao  - f;  )* 

j“l 


■ 


• • f f . (1) 


But  i). 


*r^^r  ^r  ^ ^ ^r) 

+ ap_j^(aj,  Sgpwj^  + ao  Sjwj) 

■♦■  ao(Up  Sp  '*'•.•••  +3o) 

- 2ap  2 fj  - gajw]^  2 fj  j - . ; . . {-2ao  2fj 

fj: 


♦ 2(fi)’ 


' ap-2ap  + 2(fj) 

So  that  “ j^p  •.  Sp# 


But.  from  the  conditions  (1)  above,  it  follows  that 


far  » 

%r-2  •*•••'■  ^p-1 

Sgp  • • > • . Sp 

- 

S * * * * . So 

Sp  • * • • So 

.13. 


3(n+r-l)  3(n-i^-l)[3(  [3(n-l)]*  5(8r+l) 


[3(n-l)]^  3(2r-l)  3(n+r)3(n-r-2)[3(  r) 

(n-r-1).',  (2r)!  (2r+l)J 
(n+r).'  {r.',)"* 

r 

Again, 'the  qu.anti'ty  Ai  for  the  f's  as  defined,  is 


j = r+l 

2r.'(2r+l)!(n-i“l).'  or  j-n-r 

2 

(n+r)J  (r.’)®  j»l 


(n-j).'. 


(n-j-r)l  (j-1).'  ( r-.j  >1).' 


and  the  smnination  is  in  all  cases  ufiity.. 


(n-r-1). '(2r)J.(2r+l)i 

" C - — I 

ta*r  (n+r)I  (r.-,)* 


2rl(2r+l).'(n-r-l) 
(n+r).'.  ( rJ )® 


(n+r)J  r.',“ 
(2r+l).' (n-r-1).';  2rJ 


Q.  E»  D. 


p.»f  5.  I { , 


Let  m+j  = k and  n+r  = s 
Then  L.H.  S.  = 


We  can  su;n  at  once  over  all  j , for 


f (.^)(.-o)  ■ (.t.) 


.14. 


Now  ( 2 
Ir+k- 


r \/ k-1  \ / l-l\ 

■-lj[  r J[  r J 


2r)(  k-1  V 1-1 

r Ar+k-Z /\r+l-k 


At  this  stage  it  is  usefuj.  to  pu,t  *k=Z+d 
L.H.S. . f;) . I 

Carrying  ou,t  the  suinraation  first  over  all  d,  ■ we  get 
which  may  be  written  as 

Now  it  may  be  shown  that  2 (-1)*^  (^r^)('^'^)(r+d)  ' 
for  all  a > r,  b > r. 

Putting  a * l-l,  b«s-Z,'  these  ihe<ii,alitijes  are  seen  to  apply.' 
Consequiently  we  are  left  wilth 


f 


Mf  r‘ )('"){>} 


and  this  reduces  at  once  to 


s "j 
8r+l  J 


which  = R.H.S. 


[NOTE.  We  have  omitted  the  limits  of  su(mnatibn,  since  the  summation 
is  required  for  all  valu.es  of ‘the  parameters  giving  nonzero  contribution^  ■ 
and  may  in  fact  be  taken'  from  -®  to  for  each  parameter]. 


• 15. 


TABLES  I and  lA.; 


Weilghtihg  Factors  to  be  Applied  to  Aj  in  Calcu^.atihg  ttlp'  and 
Su[ns  of  Weighting  Factors. 


Functional  Valujes  given  at  3 points 


w.  f. 


1 

1 

1 


llll  ^2 

^ 1 
2 


SUMS 


Functional  Valu.es  given  at  4 poihts 


w.  f.  for  ttlo  till  If'2 


1 

1 

1 

1 


SUMS 


3 

4 
3 

10 


3 

3 


Functional  Valu.es  given  at  5 poihts 


w.  f . fo  r lAo 


till  Iri2 


in*- 


1 

1 

1 

1 

1 


4 

6 

6 

4 


6 

9 

6 


SUMS 


20  21 


Functional  Valu.es  given'  at  6 poihts 


V.  f.  for  tHo  till  lHa 


1 

1 

1 

1 

1 

1 


5 

8 

9 

8 

5 


10 

18 

18 

10 


Ilia 


10 

10 

10 


5 

5 


in>B 


SUMS 


35  56  36  10 


TABLES  I i I A Contihujed 


Punctibnal  Valu.es  given  at  7 points 


w.  f.  for  BVo  nil 


1 

1 

1 

1 

1 

1 

1 


6 

10 

12 

12 

io 

8 


SUMS 


nil 

hi. 

m.'. 

hi. 

in. 

15 

20 

30 

15 

26 

40 

25 

6 

1 

30 

40 

15 

6 

15 

20 

126 

120 

55 

12 

1 

Rinctibnal  Valujes  given  at  9 points 


w. 'f.  for  IIVo  till  Ills 


1 

1 

1 

1 

1 

1 

1 

1 


7 

12 

15 

16 
15 
IS 

7 


21 

45 

60 

60 

45 

21 


35 

80 

100 

30 

35 


K- 


35 

75 

75 

35 


ItW 


SUMS 


8 4 252  330  220 


21 

36 

21 


78 


7 

7 


14  * 


Functibnal  Valu.es  given  at  9 points.' 


w.  f. for  tiio  nil  nil 


1 

1 

1 

1 

1 

1 

1 

1 

1 


8 

14 

18 

20 

feo 

18 

14 

8 


23 

63 

90 

100 

90 

63 

28 


iRa 


DO 

140 

200 

‘200 

140 

56 


ni*  nis  HI0  HI7  ■ 


70 

175 

225 

175 

70 


56 

126 

126 

93 


28 

49 

28 


SUMS 


9 120  462  792  715  364  105 


16 


17. 


TABLES  I (?  lA  Contihu.ed 


Functional  Valu.es  given  at  10  points 


w.  f . fo  r 


lUo 

In, 

IU2 

III3 

III* 

lUs 

Hie 

III7 

llle 

1 

, 9 

1 

36 

84 

•126 

18 

1 

84 

21 

£24 

350 

126 

34 

1 

24 

126 

1 

150 

350 

525 

336 

196 

36 

9 

25 

150 

400 

525 

441 

196 

64 

9 

1 

24 

1 

128 

350 

350 

336 

64 

36 

1 

21 

34 

224 

128 

126 

16 

36 

34 

1 

9 

1 

10 

165 

792 

1716 

2002 

1365 

560 

136 

18 

Functional  Vcilues  given  at  11  points 


w.  f.  for  tHo  tHi  rtlj 


tll4- 


ttl,  tile 


45 
108 

168  ^ Dou 

210  1050 

^5  1225 

210  1050 

168  630 

lOB  ^ 210  252 

45 


210 

630  252 

rrc/s 


210 

588 

784 


120 


f un 

588 

210  ^20 


45 

91 

45 


10 

10 


SUMS  11  220  1287  3432  5005  4368  2380  " 316  171 


20 


Fupctibnal  Valu.es  givai  at  12  points 
w.  f.  for  tllo  ttli  Ills  tile  til*'  Ills 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


11 

20 

27 

32 

35 

36 
85 
32 
27 
20 
11 


55 

135 

216 

280 

315 

315 

230 

216 

135 

55 


420 

.alb 

1120 

i^ib 

1120 

84C 

480 

165 


330 

1050 

1990 

2450 

2450 

1890 

1050 

330 


462 

1512 

2646 

3136 

2646 

1512 

482 


462 

1470 

2352 

2352 

1470 

402 


330 

960 

1296 

960 

330 


III7  In®  llle  ■ 


165 

405 

405 

165 


55 

100 

55 


SUMS  12  236  2002  e-lSS  11440  12376  3568  3676  1140  21Q- 

. 18. 


11 

11 


22 


TABLES  L iS  lA  Oontihu.ed 


Functional  Valu.es  given  at  13  points 


w. f. for  ttlo  Itli 


ttl*- 


in. 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


12 

22 

SO 

36 

40 

43 

42 

40 

S6 

30 

82 

12 


66 

.165 

270 

360 

420 

441 

420 

360 

270 

165 

66 


820 

650 

1200 

1680 

1960 

1960 

1680 

1200 

630 

220 


495 

1650 

3150 

4410 

4900 

4410 

3150 

1650 

495 


' 792 
2772 
5292 
7056 
7058 
5292 
2772 
792 


924 

3234 

5880 

7056 

5880 

3234 

924 


792 

2340 

4320 

4320 

2640 

792 


495 

1485 

2025 

1435 

495 


220 

550 

550 

220 


SUMS  13  364 


3003  11440  24310 


31824  27132  15504  5985  1540 


66 

121 

66 


253 


Functional  Valu.es  given  at  14  points 


w.  f.  fo  r ftlo 

tn.1 

riVj 

lUg 

ttl*' 

IHg 

1^7 

lUe- 

ITI9 

FVio 

1 

13 

78 

1 

286 

715 

24 

198 

1 

1 

83 

330 

880 

2475 

1287 

1716 

1716 

40 

1650 

^50 

4752 

6468 

1287 

1 

450 

9702 

6336 

715 

^6 

45 

2400 

12936 

4455 

1 

540 

7350 

11880 

48 

2940 

14112 

17640 

7425 

2200 

726 

1 

583 

8820 

15876 

14400 

49 

3136 

17640 

7425 

3025 

726 

1 

588 

8020 

14112 

31380 

48 

2940 

12036 

4455 

2200 

286 

1 

540 

2400 

7350 

9702 

63K 

45 

6468 

1237 

715 

1 

40 

450 

165d 

4950 

4752 

1716 

1 

330 

2475 

1716 

33 

880 

715 

1287 

1 

190 

24 

78 

236 

1 

1 

13 

SUMS  14 

455 

4368 

19448 

48620 

75592 

77520 

54264 

26334 

8855 

2024 

19.; 


TABLES  I 4 I A ' Contihupd 


w,  f.  for 


SUMS 


w.  f.  for 


SUMS 


Functional  Valu.es  given  at  15  points 


)Tlo  tilt  tlVa  Ills 


to*- 


tile  IHb  III7 


ttle 


ttl. 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


14 

S6 

86 

44 

50 

54 

55 
.'56 
54 
50 
44 
86 
26 
14 


91 

234 

806 

550 

675 

756 

784 

756 

675 

550 

396 

234 

91 


tS4 

1144 

2200 

3300 

4200 

4704 

4704 

■4200 

3300 

2200 

1144 

864 


1001 

3575 

7425 

11550 

14700 

15876 

14700 

11550 

7425 

3575 

1001 


2002 

7722 

16632 

25372 

31752 

3i752 

25872 

16632 

7722 

2002 


3003 

12012 

25872 

38803 

44100 

39308 

25872 

12012 

3003 


. 8482 
13723 
28512 
^600 
-39600 
28512 
18729 
3432 


3003 

11533 

22275 

27225 

22275 

11583 

3003 


2002 

.7150 

12100 

1210C- 

7150 

2002 


15  560  6188  31824  92378  167960  203400  170544  100947  42504 


Functional  Valujss  given'  at  16  points 


tAo  till  tlla  tllii 


tA*'  lAe 


tAe  lAa 


lAe  tAe 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


15 

28 

39 

48 

55 

60 

63 

64 
23 
30 
55 
48 
39 
23 
15 


105 

273 

468 

660 

8-25 

945 

1008 

1009 

945 

825 

660 

488 

273 

105 


455 

1456 

2860 

4400 

5775 

6720 

7056 

6720 

5775 

4400 

2860 

1456 

455 


1365 

5005 

10725 

17325 

23100 

26460 

26460 

231D0 

17325 

10725 

5005 

1305 


3003 

12012 

27027 

44S52 

58212 

63504 

58212' 

44352 

27027 

12012 

3003 


5005 

21021 

48048 

77616 

97020 

97020 

77616 

48048 

21021 

5005 


6435 

27456 

61770 

95040 

109900 

95040 

81776 

27456 

6435 


6435 

27027 

57915 

81675 

81675 

57915 

27027 

6435 


5005 

20020 

50325 

■©400 

39325 

20020 

5005 


16  680  8568  50388  167960  352716  407420  400314  8461D4  177100 


.20.; 


1001 

3146 

4356 

3146 

1001 


12650 


tAio 


3003 

11011 

18878 

18376 

11011 

3003 


65780 


TABLES  I & lA  Continued 


Functional  Values  given  at  17  points 


lUo 

for 


III*  Itls 


IHe 


Hi, 


Hie  Hie  Hljo 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


16 

30 

42 

52 

60 

66 

70 

72 

72 

70 

66 

60 

'52 

42 

30 

16 


120 

315 

546 

780 

990 

1155 

1260 

1296 

1260 

1155 

990 

780 

546 

315 

120 


560 

1820 

3640 

5720 

7700 

9240 

10080 

10080 

9240 

7700 

5720 

3640 

1820 

560 


1820 

6825 

15015 

25025 

34650 

41580 

44100 

41580 

34650 

25025 

15015 

6825 

1820 


4368 

18018 

42042 

72072 

99792 

116424 

116424 

99798 

72072 

42042 

18018 

4368 


8008 

35035 

84084 

144144 

194040 

213444 

194040 

144144 

84084 

35035 

8008 


11440 

51480 

123552 

805920 

261360 

261360 

205980 

123552 

51480 

11440 


12870 

57915 

135135 

212355 

245085 

212356 

135185 

57915 

12870 


11440 

50050 

110110 

157300 

157300 

110110 

50050 

11440 


8008 

33033 

66066 

81796 

66066 

33033 

8008 


SUMS  17  816  11628  77520  293930  705432  1144066  1307504  1081575  657800  296010 


Functional  Values  given  at  18  points 


nio  nil 
for 


Hie  Hie  HI4 


HI5  HIb 


Hi, 


Hie  Hie  Hlfo 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


17 

32 

45 

56 

65 

72 

77 

80 

81 

80 

77 

72 

05 

56 

45 

32 

17 


136 
. 360 
630 
910 
1170 
1386 

1540 

1620 

1620 

1540 

1386 

1170 

910 

630 

360 

136 


680 

2240 

4550 

7280 

10010 

12320 

13860 

14400 

13860 

12320 

10010 

7280 

4550 

2240 

680 


2380 

9100 

20475 

35035 

50050 

62370 

69800 

69300 

62370 

50050 

35035 

20475 

9100 

2380 


6188 

26208 

63063 

112112 

162162 

199584 

213444 

199584 

162162 

112112 

63063 

26208 

6188 


12376 

56056 

140140 

252852 

360360 

426888 

426888 

360860 

252252 

140140 

66C66 

12376 


19448 

91520 

231660 

411840 

566280 

627264 

566280 

411840 

231660 

91520 

19448 


24310 

115830 

289575 

495495 

637065 

637065 

495495 

289575 

115830 

24310 


24310 

114400 

275275 

440440 

511225 

440440 

875875 

114400 

24310 


19448 

88088 

198198 

286286 

286286 

198198 

88088 

19448 


SUl-lS  18  969  15504  116280  497420  1352078  8496114  S26R760  S124550  2220075  1184C40 
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TABLES  I & I A Continued 


Fimctional  Values  given  at  19  points 


w.  f. 

for  Itlo  lUl  Itl2  tlia  ^4 


1U7 


trie 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


18 

34 

48 

60 

70 

78 

84 

88 

90 

90 

88 

84 

78 

70 

60 

48 

34 

18 


153 

408 

720 

1050 

1365 

1638 

1848 

1980 

2025 

1980 

1848 

1638 

1365 

1050 

720 

408 

153 


816 

2720 

5600 

9100 

12740 

16016 

18480 

19800 

19800 

18480 

16016 

12740 

9100 

5600 

2720 

816 


3060 

11900 

27300 

47775 

70070 

90090 

103950 

108900 

103950 

90090 

70070 

47775 

27300 

11900 

3060 


8568 

37128 

91728 

168168 

252252 

324324 

365904 

365904 

824324 

252252 

168168 

91728 

37128 

8568 


18564 

86632 

224224 

420420 

630630 

792792 

853776 

792792 

630630 

420420 

224224 

86632 

18564 


31824 

155584 

411840 

772200 

1132560 

1859072 

1359072 

1132560 

772200 

411840 

155584 

31824 


43758 

218790 

579150 

1061775 

1486485 

1656369 

1486485 

1061775 

579150 

218790 

43758 


48620 

243100 

629200 

1101100 

1431430 

1431430 

1101100 

629200 

243100 

48620 


43758 

213928 

528528 

858858 

1002001 

858S58 

528528 

213928 

43758 


SUMS  19  1140  20349  170544  817190  2496144  5200300  7726160  8436285  6906900  4292145 


Functional  Values  given  at  20  points 


w.  f. 

for  ttlc  till  titj  tlia  tH*  Me  Me 


M7  Me 


Me  Mi'o 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


19 

36 

51 

64 

75 

84 

91 

96 

99 

100 

99 

96 

91 

84 

75 

64 

51 

36 

19 


171 

459 

816 

1200 

1575 

1911 

2184 

2376 

2475 

2475 

2376 

2184 

1911 

1575 

1200 

816 

450 

171 


909 

3264 

6800 

11200 

15925 

20384 

24024 

26400 

27225 

26400 

24024 

20384 

15925 

11200 

6800 

3264 

969 


3876 

15300 

35700 

63700 

95550 

126126 

150150 

163350 

163350 

150150 

126120 

95550 

03700 

35700 

15300 

3876 


11628 

51408 

129948 

244808 

378378 

504504 

594594 

627204 

504594 

504504 

378378 

244608 

129948 

51408 

11623 


27132 

129948 

346528 

672672 

1051050 

1387386 

1585^4 

1585584 

1387386 

1051050 

072672 

346588 

129948 

27138 


50388 

254592 

700128 

1372800 

2123550 

2718144 

2944656 

2718144 

2123550 

1372800 

700128 

254502 

50388 


75582 

393822 

109^50: 

21.23550. 

3185325 

3864861 

3864861 

3185325 

2123550 

1093950 

393822 

75582 


92378 
486200 
1337050 
2516800 
3578575 
4008004 
3578575 
2516800 
1337050 
486200 
92378 


92378 

481338 

1283568 

2290^8 

3006003 

3D06003 

2290288 

1283568 

481333 

92378 


SUMS20  1330  26334  245157  13C7504  4457400  10400600  17383860  21474180  20030010  14307150 
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TABLE  II 


Summary  Table  of  Sums  of  the  Weighting  Factors. 

S 
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CO 
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3 

CO 

1 

1 

General  form:  The 

sum  of  the  weighting  factors 

2 1 

for  Itlj,,  when  functional  values  are  given  at  n 

3 4 

1 

points,  is 

the  binomial  coefficient  1 

n+r  1 

I • 

4 10 

6 

1 

\ 

,2r+iy 

5 20 

21 

8 

1 

6 35 

56 

36 

10 

1 

7 56 

126 

120 

55 

12 

1 

8 84 

252 

330 

220 

78 

14 

1 

9 120 

462 

792 

715 

364 

105 

16 

1 

10  165 

792 

1716 

2002 

1365 

560 

136 

18 

1 

11  220 

1287 

3432 

5005 

4368 

2380 

816 

171 

20 

1 

12  286 

2002 

6435 

11440 

12876 

8568 

3876 

1140 

210 

22 

13  364 

3003 

11440 

24310 

31824 

27132 

15504 

5985 

1540 

253 

14  455 

4368 

19448 

48620 

75582 

77520 

54264 

26334 

8855 

2024 

15  560 

6188 

31824 

92378 

167960 

203490 

170544 

100947 

42504 

12650 

16  680 

8568 

50388 

167960 

352716 

497420 

490314 

346104 

177100 

65780 

17  ■ 816  11628 

77520 

298930 

705432 

1144066 

13C7504 

1081575 

657800 

2^010 

18  969  15594  116280 

497420 

1352078 

2496144 

3268760 

3124550 

2220075 

1184040 

19  1140  20349  170544 

817190  2496144 

5200300 

7726160 

8^-36285 

6906900 

4292145 

20  1330  26334  245157 

1307504  4457'100  10400600  17383860  21474180  20030010 

14307150 

TABLE  III 


- (?>j.  is  given  by  Kltl*,  where  the  constant  K is  equal  to  the  sum 
of  the  weighting  factors  for  Itlj.,  given  in  the  table  above,  divided  by  the 

binomial  coefficient  ^ function  of  r only  and  independent  of  n. 

This  is  tabulated  below  for  r up  to  10, 

r01234567  8 9 10 

1 2 8 20  70  252  924  3432  12870  49620  184756 
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